Negative group delay for Dirac particles travelling through a potential well 
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The properties of group delay for Dirac particles travelling through a potential well are investi- 
gated. A necessary condition is put forward for the group delay to be negative. It is shown that this 
negative group delay is closely related to its anomalous dependence on the width of the potential 
well. In order to demonstrate the validity of stationary-phase approach, numerical simulations are 
made for Gaussian-shaped temporal wave packets. A restriction to the potential-well's width is 
, obtained that is necessary for the wave packet to remain distortionless in the travelling. Numerical 

comparison shows that the relativistic group delay is larger than its corresponding non-relativistic 
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I. INTRODUCTION 
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The question of how much time it takes quantum particles to tunnel through a potential barrier has been controver- 
sial for these decades 0, 0, 01 ■ Theoretical investigations 0, IS IE S 13 an d experimental researches 0, El 0, 0, H^ ? 
] show that the group delay for some kinds of barriers, also known as the phase time in the literature, which describes 
the motion of a wave pack et p eak |l4j . has the well-known superluminality. In addition, the faster-than-light propa- 
gation was also predicted |15| | and experimentally verified |16l H71 Il8| for light pulses through anomalous dispersion 
' media. In a previous paper |19| . Li and Wang have elaborated the superluminal and even negative properties of 
Ck the group delay for quantum particles travelling through a potential well, instead of tunnelling through a potential 
+-> ' barrier. This counterintuitive phenomenon resulted from the interference of multi-reflected waves in the potential well 
was demonstrated in a microwave analogy experiment |2f| . Recently, the concept of negative group delay has been 
extended to microelectronics plj . 
q-( Most of the theoretical works on tunnelling times in quantum mechanics rely on Schrodinger's non-relativistic 
theory. Such a theory has apotential deficiency in accurately addressing the question of causality [2^, [2j| . For this 
reason, a few authors |22, J23, 0, 113 extended the analysis of the tunnelling to Dirac's fully relativistic quantum 
theory. Leavens and Aers [24J used the stationary-state method to analyze the Larmor-clock transmission times for 
single barriers and resonant double-barriers. Krekora et al. p2| solved numerically the time-dependent Dirac equation 
for a quantum wave packet tunnelling through a potential barrier. And Petrillo and Janner |25j studied the dynamics 
of wave-packet tunnelling through a barrier. In a recent work j23j |. we discussed an energy-transfer associated traversal 
time for Dirac particles tunnelling through a potential barrier. 

The purpose of this paper is to investigate the properties of the group delay for Dirac particles travelling through 
a potential well. It is shown that it behaves superluminal and even negative in much the same way as in the non- 
relativistic quantum mechanics |19| . The negativity of the group delay is closely related to its anomalous dependence 
on the width of the potential well around transmission resonances. In order to demonstrate the validity of the 
stationary-phase approximation, numerical simulations are made for Gaussian-shaped temporal wave packets. A 
restriction to the width of the potential well is given that is necessary for the wave packet to remain distortionless 
in the travelling. Finally, a numerical comparison shows that the Dirac's relativistic group delay is larger than its 
corresponding non-relativistic one. 
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II. RELATIVISTIC GROUP DELAY AND ITS NON-RELATIVISTIC LIMIT 



Consider Dirac particles of precisely defined incident energy E and of helicity +1, travelling through a one- 
dimensional rectangular potential well — Vo<d(z)Q(a — z) (with Vq positive, representing the depth of the potential 
well). Let the incident wave function be 



4>in(z) = 



( 1 \ 



hkc 

V o j 



e ikz ,(z<0), 



(1) 



where k = (E 2 — juV) 1 / 2 /he, fi is the mass of incident particles, and c is the speed of light in vacuum. Then Dirac 
equation and boundary conditions give for the transmitted wave function 



ip tr (z) = F 



( 1 \ 



hkc 

V o j 



e lk(z - a \(z > a), 



(2) 



where the transmission coefficient F — e 1 ^ / f is determined by the following complex number, 

fe** = cosk'a + (i/2)( X + l/x) sinfc'a, 

so that 



/ k a 1 i 
= mt I 1 — ) 7r + tan 

7T 2 



1 



f 



- ( x H — ) tanfc'a 



X 



(3) 



int(-) stands for the integer part of involved number, k' = [(E + Vb) 2 — (Pc 4 } 1 / 2 /he. Noting that the real parameter 
X is defined as 



X 



k E + V + fic 2 
k' E + fic 2 



and has the property that < x < 1 ■ Obviously, </> is the phase shift of transmitted wave J2J at z — a with respect to 
the incident wave Q at z = 0. The transmission probability T is a periodical function of the width a of the potential 
well, 



T 



4 X 2 



p 4^2 + ^2 _ 1)2 gm 2 k , a 



(4) 



by 



The group delay, defined as the derivative of the phase shift <fi with respect to particle's energy E 0, ^| , is given 



T4, 



T 



dE 2 x hk'c 



H + X>){E + Vo)-{l-xY M2E + V ° )iihl2kla 



h 2 k 2 



2k'a 



(5) 



For the following comparisons, let us first give its non-relativistic limit. By non-relativistic limit it is meant that the 
speed of incident particles is much smaller than c, so that their kinetic energy E' = E — fic 2 is much smaller than 
their rest energy fic 2 , E' -C /ic 2 . It is also meant that the interaction energy Vb satisfies Vq <C ■ In this limit, we 
have k y / 2~JlE 7 /h, k' w [2(J,(E' + V )] 1/2 /h, and \ ~ k/k'. Collecting all these, we get 



T 4> ~ T 



2/ia k 2 (k 2 + k' 2 )/k 4 - sin2k'a/2k'a 



hk Ak 2 k' 2 /k 4 + sm 2 k'a 

which is expected from Schrodinger's non-relativistic theory [Tsj. where fcg w ^/2pVo/h. 



(6) 
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III. NEGATIVE PROPERTY OF THE GROUP DELAY 



In this section, we discuss the negative property of the group delay. It is seen from Eq. J5J that when inequality 

(l + X 2 )(E + V )<(l-x 2 ) 1 - 



j - , , , 2^ fjV (2E + V ) sin 2k' a 
' Wk? 2k'a 



holds, the group delay is negative, < 0. Since sm2k'a/2k'a < 1, the above inequality leads to the following 
necessary condition for the group delay to be negative, 



(l + X A )(E + V Q ) < (l-xl 



2 .HV (2E + V ) 



h 2 k 2 



(7) 



which can be expressed as a restriction to the total energy of incident particles as follows, 



E<E t 




1/2' 



1/3 



1/2' 



1/3 



(8) 



This means that when the energy of incident particles satisfies Eq. (JSJ) , that is to say, the energy of incident particles 
E is less than a threshold energy E t , one can always find a width a of the potential well at which the group delay 
is negative. Of course, Et in Eq. JHJ is always real and larger than /ic 2 . In the case of Vb < 2/ic 2 /3-\/3, E t can be 
rewritten as 



E t = /ic 2 



Vo 
2/icr 



+ i 



Vo 
2/^c 2 



1/2' 



1/3 



C.C. 



(9) 



In the non-relativistic limit, Vb -C /^c 2 , we obtain E t ~ /ic 2 + Vn/2. Therefore, the necessary condition © now reduces 
to E' < Vq/2, as observed previously in Schrodinger's theory 

Fig. Q]shows a typical example of the dependence of on the width a, where the well depth Vq = OAfic 2 > 2/xc 2 /3v / 3 
(E t = 1.16^ic 2 ), the total energy E = l.Ol^c 2 < E t , and a is re-scaled to be k'a. For comparison, Fig. ^also shows 
the periodical dependence of transmission probability T on a under the same conditions. It is interesting to note 
that the oscillation of the group delay with respect to a is closely related to the periodical occurrence of transmission 
resonances at k'a = mm (m = 1,2,3...). 
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FIG. 1: Dependence of the group delay (in the unit of h/^c 2 ) on potential-well's width a, where E — 1.01/ic 2 , Vo = 0.4/ic 2 , and 
a is re-scaled to be k'a. Here the periodical dependence of transmission probability T on a is also depicted by dashed curve. 



On the one hand, at resonances, the group delay becomes 

1 



1\ E + V a a 
hk'c c c ' 



(10) 
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which is proportional to a. Its corresponding group velocity is less than the velocity of light in vacuum, c. On the 
other hand, the derivative of group delay ta, with respect to a is, at resonances, 



<9t0 
da 



1 



k' a—mir 



2 X hk'c 2 



(l + X 2 )(E + V )-(l-x 2 ) 



fiV Q (2E + V a ) 



h 2 k 2 



(11) 



When the necessary condition JJJ is satisfied, it is negative. This shows that the group delay depends anomalously 
on the width around resonance points. In other words, it decreases with increasing the width of the potential well at 
resonances, as displayed clearly in Fig. ^ 

In connection with the negative characteristics of the group delay, the phase shift J3J shows a quantum-like behavior 
with respect to the potential-well's width. Fig. [21 indicates such a behavior, where E = 1.01/xc 2 and Vo = 0.4/xc 2 . 
This quantum-like behavior is also related to the periodical occurrence of the transmission resonances. It can be seen 
from Eq. © and Fig. [21 that at resonances, k'a = rn7r, the phase shift becomes cj) = k'a and changes rapidly around 
here with respect to a. However, when the width is far from the resonance points as is in the middle between two 
adjacent resonance points, the phase shift changes slowly. 




k 1 a 

FIG. 2: The quantum-like behavior of phase shift © with respect to the width of potential well, a, where E = 1.01/iC 2 and 
Vo = 0.4/jc 2 



In addition, it is also indicated from Eq. JSJ that the group delay depends not only on the potential-well's width a, 
but also on the incident energy E and the depth of the potential well Vq. To see the latter more clearly, we convert 
it to a dimcnsionless form. Denoting a — E//j,c 2 , (3 — Vo/fic 2 , and 7 = h/afic, we have 



TO 



T 



k'a 



2x {a + PY 



1 



(3(2a + f3) sin2fc'a 
a 2 -I 2k'a 



(12) 



where tq — h/fic 2 , k'a — [(a + 0) 2 — l] 1//2 /7- When the kinetic energy of incident particles E' is small enough that 
E' I /xc 2 — + with the depth of the potential well remaining finite, we have a — * 1, so that x ~ * 0- I n this limit, the 
group delay (|T2l has the following form, 



lim T ± = 

<*-»! tq 



' 13 + 2 

{a 2 - 1)13 



cot k'c 



which approaches negative infinite when cot k'a > 0. Fig. shows such a dependence of the group delay on the 
incident energy E, where (3 = 0.4 and 7 = 0.01. A strange phenomenon occurs here that for a given potential well, 
the absolute value of the negative group delay becomes larger when decreasing the incident kinetic energy. Of course, 
the transmission probability in this limit tends to zero in the following way, 



lim T 



4x 2 



4% 2 + sin k'a 



(13) 



so that very few particles can travel through the potential well at this negative group velocity. 

In order to demonstrate the validity of the above stationary-phase method in this problem, we proceed to numerical 
simulations of the group delay for a Gaussian-shaped wave packet. The incident wave function of Dirac particles is 
assumed to be, at z = 0, 



*m(*)U=0 = 



1 



Hkpc 





\ 



exp(-< 2 /2w 2 - iE t/H), 



(14) 
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FIG. 3: Dependence of the group delay (in the unit of h/fic 2 ) on the incident energy a, where @ = 0.4 and 7 = 0.01. 



which has the Fourier integral of the following form 

*in(*)U=0 ' 



'2tt 



A(E)ip(E ) exp(-iEt/h)dE, 



(15) 



where k = {El - /^c 4 ) 1 / 2 /he, the spinor %P(E ) is [1,0, hk c/{E + fic 2 ),0] T , the energy spectral distribution A(E) 
with the central energy Eq is given by 

A{E) = [w/h) ex P [-(w 2 /2h 2 )(E - E ) 2 }, 

and w is the temporal width of the wave packet l)14|l . The transmitted wave function takes the following form, 

1 



/2tt 



F(E)A(E)ip{Eo)exp{i[k(z - a) - Et/h]}dE, 



The numerically calculated group delay, tv , is defined here by 



% r (z = a,T*)\ 2 = max{|4< tr (z = a,t)\ 2 }. 



(16) 



(17) 



Since the incident wave is of perfect Gaussian shape, the integral limit in Eq. (|15J) and hence in Eq. Ijl6|l should be 
from —00 to +00. But the energy of incident particles must be larger than /ic 2 . So the real integral in numerical 
simulations is taken to be from \m? to +00. 

Calculations show that the stationary-phase approximation (JEj) for the group delay is in good agreement with the 
numerical result, especially when the energy spectral distribution is sharp. In Fig. 0] we show such a comparison 
between theoretical and numerical results, where Eq = 1.01/ic 2 , Vq = 0.4/ic 2 , and the temporal width w = 300ro. For 
the chosen temporal width, the corresponding energy spreading AE = h/2w = fie 2 /600 is narrow enough that the 
integral in Eq. I|16fl can be performed from \ic 2 to 2/ic 2 without changing the result significantly. 

As pointed out in the optical analog |26| . for an incident wave packet of energy spreading AE, the corresponding 
spreading of k'a should be much smaller than tt, the period of \F\, in order that the stationary-phase approximation 
is valid. With the energy spreading AE = h/2w, this leads to the following restriction to the width of the potential 
well, 



a <C 2ttw 



dE 
hdk' 



(18) 



It is noted that dE/hdk' is nothing but the group velocity of particles in the region of potential well. By introducing 
a characteristic length L which is defined as wdE/hdk', the above restriction is simplified to be a <C 2irL. With this 
restriction, the temporal wave packet can travel through the potential well with negligible distortion. 

All the above discussions show that the group delay (JSJ) in Dirac's relativistic theory can be superluminal and even 
negative in much the same way as in Schrodinger's non-relativistic theory |19| . In this sense, the superluminal and 
negative properties is not an artifact due to the deficiency of non-relativistic quantum theory. Rather, it is a real 
effect. In the next section, we investigate the difference between the group delays in relativistic and non-relativistic 
theories. 
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FIG. 4: Comparison of theoretical and numerical results of group delays (in the unit of h/^c 2 ) with respect to a, where 
E Q = l.Ol^ic 2 , Vo = OAfic 2 , w = 300r and a is re-scaled to k'a. The theoretical result for the group delay 10 is depicted by 
the real curve, and the dash curve corresponds to the numerical result for the group delay defined by 11711 . 



IV. RELATIVISTIC EFFECT ON THE GROUP DELAY 



As mentioned in Sec. II, when the incident kinetic energy E' <§; [m? and the potential- well's depth Vo <C ^c 2 , the 
group delay tends to its non-relativistic limit JJ]J. To show their differences, we rewrite the non-relativistic group 
delay in the following dimensionless form, 



k'a 



(a - l)(2a + P - 2) - (3 2 sin2k'a/2k'a 



l)(a + /3- 1) 4(a -l)(a + (3- 1) + /? 2 sin 2 k'a 



(19) 



where a, /3, and 7 are the same as before, k'a = [2 (a + f3 - l)] 1/2 /7 in the non-relativistic quantum mechanics. 

In order to address the relativistic effect on the group delay, we draw in Fig. \E\ the dependence of relativistic and 
non-relativistic group delays on the width of the potential well for a = 1.01 and (3 — 0.2 (Vo = 0.2/ic 2 < 2/ic 2 /3v / 2), 
where the relativistic group delay is shown by solid curve, the corresponding non-relativistic one is shown by dashed 
curve, and the width of potential well is re-scaled by their own k' to be k'a. 




FIG. 5: Comparison of relativistic and non-relativistic group delays (in the unit of ft//ic 2 ) for a = 1.01 and f3 = 0.2. The real 
curve denotes the relativistic group delay H12I . and the dashed curve corresponds to the non-relativistic group delay 11911 . 



From Fig. we see that the relativistic group delay for travelling through a potential well is larger than the 
corresponding non-relativistic one. This is in agreement with the result of Leavens and Aers [24|. who observed 
that the local mean velocity for transmitted particles is reduced due to relativistic effect for resonant double barriers 
where the group delay is always larger than zero. The difference between the two group delays are very small for the 
low energy and potential-well depth. It is expected that the agreement between the group delays in relativistic and 
non-relativistic theory is obtained when E' <C /xc 2 and Vo <C /ic 2 . 
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V. CONCLUSIONS 



In summary, we have investigated the negative property of the group delay for Dirac particles travelling through a 
quantum potential well. A necessary condition Q is given for the group delay to be negative, which is a restriction 
on the energy of incident particles and reduces to E' < Vq/2 in non-relativistic limit. The relation of the negativity of 
the group delay with its anomalous dependence on the width of the potential well around transmission resonances is 
discussed. In order to demonstrate the validity of the stationary-phase approach, numerical simulations are made for 
a Gaussian-shaped temporal wave packet. A restriction to the width of the potential well is given that is necessary for 
the wave packet to remain distortionless in the travelling. Comparison of the relativistic and non-relativistic group 
delays is also made. It is found that the value of the relativistic group delay is larger than that of the non-relativistic 
one. It should be pointed out that the negative group delay discussed here is not at odds with the principle of causality. 
When a wave packet travels through a potential well, the boundary effect leads to a reshaping of the edge of the packet 
|llll27j | which results in an effective acceleration of the wave packet. This phenomenon is similar to but different from 
the negative group-velocity propagation performed by Wang et. al. 18] where the gain-assisted anomalous dispersion 
plays the role. We hope that this work may stimulate further experimental researches in electronic domains. 
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